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Abstract 

We consider the weak localization correction to the conductance of a ring connected 
to a network. We analyze the harmonics content of the APtshuler-Aronov-Spivak (AAS) 
oscillations and we show that the presence of wires connected to the ring is responsible for a 
behaviour different from the one predicted by AAS. The physical origin of this behaviour is 
the anomalous diffusion of Brownian trajectories around the ring, due to the diffusion in the 
wires. We show that this problem is related to the anomalous diffusion along the skeleton 
of a comb. We study in detail the winding properties of Brownian curves around a ring 
connected to an arbitrary network. Our analysis is based on the spectral determinant and 
on the introduction of an effective perimeter probing the different time scales. A general 
expression of this length is derived for arbitrary networks. More specifically we consider the 
case of a ring connected to wires, to a square network, and to a Bethe lattice. 

PACS numbers : 73.23.-b ; 73.20.Fz ; 72.15.Rn ; 02.50.-r ; 05.40.Jc. 

1 Introduction 

At a classical level, a network made of diffusive wires can be described as an ensemble of 
classical resistances. Quantum corrections bring a small sample specific contribution whose 
disorder average is called the weak localization (WL) correction. This contribution is sensitive 
to a magnetic field and makes the average conductivity of a phase coherent ring be a periodic 
function of the magnetic flux through the ring with periodicity (/>o/2, where 4>q = h/e is the 
flux quantum. This is the so-called Al'tshuler-Aronov-Spivak (AAS) effect pQ. The case of an 
isolated ring was considered by AAS who showed that the average correction to the classical 
conductivity varies as : 

An(ft\- 6 % sinh(L/L y ) 

where L is the perimeter and 9 = Att^/^q the reduced flux. Phase breaking mechanisms are 
taken into account through the characteristic length L v called the phase coherence length. The 
harmonics of the oscillations 

r2n jn 2 

AffW = / — A<r(0) e~ ine = ~L V e -l"l L /^ (2) 

decay exponentially with the perimeter L of the ring and the order n of the harmonic 1 . The 
relation (J2J) was derived for an isolated ring and it is not clear, when studying the transport 

x Note that Act** 1 ' = Act' - ™' is a general property related to the symmetry Aa(9) = Act(— 9) and we will 
simply omit the absolute value in the following. 
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through a ring connected by wires to reservoirs, how this correction to the conductivity is 
related to the correction to the conductance. Moreover the presence of the connecting wires can 
seriously modify the harmonics content of the AAS oscillations. In this paper we show that, if 
the perimeter L of the ring is much smaller than L v and if the connecting wires are much longer 
than Lp, the nth harmonic decreases like exp —nW N a L/L,p, that is faster than (|2jl. N a is the 
number of wires attached to the ring. 

The paper is organized as follows : in section |2] we briefly recall the general expression of 
the WL correction on a network and we present some results for the conductance through a 
connected ring. The rest of the paper is devoted to the analysis of the new behaviour of the 
harmonics. This is achieved by studying in detail the winding properties of Brownian curves 
around a ring connected to a network. In section |3] we present our general method, based on the 
spectral determinant which encodes the necessary information on the network. In section 0] we 
consider specifically the question of a ring connected to one or several wires. We show that the 
new behaviour exp —n-sjNaL/L^ is the signature of an anomalous diffusion around the ring. In 
section |21 we study the case where the ring is connected to an arbitrary network. We relate the 
winding properties around the ring to the recurrent character of the Brownian motion in the 
network. 




Figure 1: A network connected to reservoirs and pierced by a magnetic flux 4>. The wavy lines 
represent connection to large contacts (external reservoirs). 



2 Quantum transport through a connected ring 

Recently, we have obtained a general expression for the weak localization (WL) correction on 
a network [2j Ej. Let us remind first how it reads for a wire. The classical conductance of a 
wire of section s and length L is given by the Ohm's law G cl = ctqs/L where o~o is the Drude 
conductivity. This result can be rewritten in terms of the total transmission through the wire, 
i.e. the dimensionless conductance, T^ ire = adN c £ e /L, where N c is the number of channels, £ e 
the elastic mean free path and a numerical constant depending on the dimension {ct\ = 2, 
Oil = tt/2 and = 4/3). The interferences between reversed trajectories are described by the 
so-called cooperon P c (x, x) which measures the return probability for a diffusive particle and 
which is solution of a diffusion equation that we will recall in section 03 These interferences give 
rise to the WL correction which is expressed as 

2 f L 

AT wire = J dxP c (x,x). (3) 

More generally, for a network attached to leads a and j3 as shown on figure ^ the classical trans- 
mission coefficient is obtained by classical Kirchhoff laws and it has the form T^a = otdN c £ e / ' C a p 
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where the equivalent length C a p is function of the lengths of the wires, the topology of the 
network and the way it is connected to external reservoirs. Note that C a p is simply proportional 
to the equivalent resistance. 

On such a network, because of the absence of translation invariance, we have shown recently 
that the cooperon must be properly weighted when integrated over the wires of the networks to 
get the WL correction to the transmission coefficient. Equation (J2J) generalizes as [2] : 

2 dT d r 

AT a/3 = — — £ -f- \ dx P c (x, x) . (4) 

The sum runs over all wires (liv) of the network. The cooperon is integrated over each wire 
(fj,f), with a weight which is simply the derivative of the classical transmission with respect to 
its length l^ u . This WL correction depends on the lengths of the wires, the phase coherence 
length L,p, the magnetic field and the topology of the network. 




Figure 2: A mesoscopic ring connected at two reservoirs. 

We consider the transport through the ring of figure |2 The classical transmission reads 
TrL* = "d^de where V~}, , = Zr 1 + Z7 1 . To simplify the calculations we consider below the 

rm g la+l c //d + k ' c//d c d 

symmetric case l a = lb and l c = id = L/2. The calculation of the WL correction AT r i ng (#), given 
by @, requires the construction of the cooperon in the network, explained in [21 01 ■ We do not 
give further details nor the full result ;3; and we only present two limiting cases for the Fourier 
decomposition AT r -"' g : 

• The weakly coherent network L„ -C L,l a . 

AT (n) ~ ( -V" e' n L/L ^ for n > (5) 

^ring- 4( 2 / a + L/4)2 lorn>u. [b) 

We get an exponential behaviour coinciding with the one predicted by AAS. The factor (2/3) 2n 
is related to the probability to cross 2n times a vertex of coordinence 3 [7J- This result shows 
that the connecting wires play an important role and that the determination of L« from ratio 
of harmonics using formula © can lead to a wrong estimate. 

• The mesoscopic ring connected to long wires L <C <C l a - 

The harmonics decay presents a strikingly different behaviour from the one given by AAS. 
This can be understood in the following way : because L <C L v , whenever a diffusive trajectory 
encircles the ring, it explores a distance in the arm which can be much longer than the perimeter 
of the ring. This leads to an effective perimeter L c g ~ ^Jlh^L 3> L which is responsible for a 
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decrease of the harmonics with n faster than (J2J). In the next sections, we examine in detail the 
effect of the arms on the winding properties around the loop in order to analyze the origin of 
the new behaviour ©. 

3 Winding of Brownian trajectories and spectral determinant 

In) 

In this section, we recall how the harmonics of the WL correction, AT ring , can be related to 
the winding properties of diffusive trajectories around the ring. For this purpose we introduce 
the probability to start from a point x and come back to it after a time t having performed n 
windings around the loop. This probability can be written with a path integral as : 

px(t)=x , 

W n (x,x;t)= Vx(T)e-y° dT± ^ 2 5 nM[x(T)] (7) 

Jx(0)=x 

where M[x(t)] = J *dri(r) is the winding number of the diffusive trajectory x(r). We have 
chosen a diffusion constant equal to unity D = 1. If we write 5 n ^ = ^e 1 ^ - "^, the 
coupling of the winding number to the conjugate parameter 9 in the action is interpreted as the 
coupling to a magnetic flux. The Laplace transform with respect to time relates the probability 
W n (x,x;t) to the cooperon P c (x,x) : 

/•oo 

/ dtW n (x,x;t)e~^ 
Jo 



The cooperon P c (x, x') is solution of a diffusion-like equation [7 — D^]P c (x, x') = 5(x — x') where 
the covariant derivative is = ^ — iA(x) with A{x) = 0/L for x inside the ring. In the 
WL theory, the cooperon describes the contribution of quantum interferences, that are limited 
by phase breaking mechanisms. In this respect, the Laplace parameter 7 plays the role of 
a phenomenological parameter that selects the diffusive trajectories for times t < I/7. This 
parameter is related to the phase coherence length L v , that gives the length scale over which 
quantum interferences can occur, by 7 = l/L^. The n-th harmonics of the WL correction AT^ 
is given by an integral over x of the Laplace transform f °° dt W n (x,x;t) e~ 7 *. The integration 
over the network is performed by weighting the wires with the coefficients given in eq. @. 

The relation between WL and properties of the Brownian motion was formulated in other 
works like [IJ[3]- 

3.1 Probability averaged over the network 

First, we do not consider the dependence on the initial point x and average over it. We define 
W n (t) = J dxW n (x, x;t). In order to study this quantity it is useful to introduce the spectral 
determinant, defined as #(7) = 1X1(7 + E n (9)), where {E n (6)} is the spectrum of the operator 
— D^. The spectral determinant is related to the cooperon in the following way : J dx P c {x, x) = 
Tr{ } = Yin 7 +g = ^ m ^(7)- Therefore the probability W n (t) can be conveniently 
written as : 

/■OO r2n jo a 

I ^)°-*=l ^-'V nS(7) - (10) 

The interest to introduce the spectral determinant is that it is a global quantity encoding all 
informations about the spectrum. This approach is especially efficient for networks thanks to 
the compact expression of 5(7) obtained in : it can be related to the determinant of a finite 



2 k 



d9 
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271 dB_ 
2^ 



e~ ind P c (x,x) . 



x(t)=x 



dte'^ / £>xe-Uo dri2+i rJo d ™ (8) 



x(0)=x 



(9) 
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size matrix that encodes the information on the topology of the networks, the lengths of the 
wires, the magnetic fluxes and the way the network is connected to external reservoirs. This 
relation is briefly recalled in the appendix 1X1 



3.2 Probability for a fixed initial condition 



It is also interesting to study the probability W n (xo, xo]t) for n windings in a time t when the 
initial condition xq is fixed. This requires some local information which is obtained from the 
construction of P c (xq,xq) inside the network (the expression can be found in jHJ). We propose 
here a method to obtain this local information still using the spectral determinant. This latter 
encodes a global information since it results from a spatial integration of the cooperon over the 
network. It is related to a sum over the eigenvalues of the Laplace operator X^n j+e • On 
the other hand P c {xq,xq) requires some local information on the eigenfunctions. In order to 
extract this local information we introduce the modified cooperon P c (Ao) , solution of [7 -D 2 X + 
Xq8(x — xq)]Pc X °\x, x') = 5{x — x'). The corresponding spectral determinant is 5 7 ln5^ A °^(7) = 



/ dxP 



(Ao), 



x,x 



Its expansion in powers of the parameter Aq leads to : <9 7 lnS , ( A °)(7) 



d 7 lnS(7) + Xod.yP c (xQ, xq) + O(Ag). The cooperon at xo can be obtained by computing the 
spectral determinant S^ ^^) for a 5 added at xq : 



P c (x ,x ) = -^lnS^\ 1 ) 



An=0 



(11) 



and inserted in (jUJ). For a network, the computation of S^ ^"/) is especially simple : Ao is inter- 
preted as the parameter involved in the boundary condition at a vertex at xq (see appendix IX)) . 



4 Winding around a ring with an arm 

In this section we study how the winding properties inside a ring are affected by the presence 
of an arm. We first recall the simple case of the isolated ring and consider the case of the ring 
with one arm (figure OK) f° r different kinds of boundary conditions at the end of the arm. 



4.1 Spectral determinant and winding properties 

A reminder : the isolated ring.— In order to illustrate our formalism, we start by recalling 
few simple results about the isolated ring, that will be useful for the following. The quantity at 
the basis of our analysis is the spectral determinant, that reads for the isolated ring : £(7) = 
2(cosh^/7"L — cos 9). The relation (jl()j) shows that the Laplace transform of the probability 

W n (t) require a Fourier transform of In 5(7) = ^ cosh^L^cos 6 ■ We obtain 

OO T _ 

dtW n (t)e-^ = -—e- n ^ L , (12) 
2 \/7 

which is the exponential behaviour recalled in eq. (f2"j). It is worth noticing that the translation 
invariance of the system implies that W n (t) = LW n (x, x; t). The inverse Laplace transform leads 
to 

I (riL) 2 

W n (x,x;t) = -^e~— (13) 



2V7rf 

which characterizes normal diffusion in the ring. We can study the scaling of the winding number 
with time by computing n t = w (n 2 )t, where (•••)* denotes averaging over all trajectories for a 
time t. The ring possesses a characteristic time scale : the Thouless time. In units such that 
D = 1, it reads tl = L 2 . It is the time needed by a diffusive trajectory to explore the ring. In 
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the short time limit t <C tl , we obtain n t ~ y/2 exp — j^- and in the long time limit t 3> rx, we 
get n t ~ \/2t 1/2 /i- 




Figure 3: Left : ^4 network with one loop and one arm. A general boundary condition (Xi) is 
chosen at vertex 1. Ai = describes a reflecting boundary (isolated network). Ai = oo descibes 
an absorbing boundary (connection to a reservoir). Right : The study of the winding in the ring 
is equivalent to study the displacement along the skeleton of a comb. 

Dirichlet boundary at the end of the arm. We now consider the winding in a ring 
with an arm. A Dirichlet boundary at the end of the wire describes absorption (connection 
to a reservoir). In terms of the parameter introduced in the appendix ^ this corresponds to 
Ai = oo. The spectral determinant of the ring of figure OK is easily obtained (see 00131 and 
appendix EJ . It presents a structure similar to the one of the isolated ring, 

^yS Dil ( 7 ) = 2 sinh ^b [cosh ^L^ T - cos 9] , (14) 

where we have introduced L^f, defined by 

cosh ^/7~L^ r = cosh ^fyL + - sinh y/jL coth ^fyb . (15) 

Interestingly, the flux dependance S(-y) oc cosh^/7~L e ff — cos 9 immediatly leads to the behaviour 
J °° dtW n (t) e~ 7 * oc e _n v / T ; - Lcff . In general the length L e g depends on 7. Since 7 is the parameter 
conjugate to the time t, it probes times t ~ I/7. Therefore L e g(7) is the effective perimeter for 
winding trajectories at a time scale t = I/7. 

In the limit 6 — > the spectral determinant becomes <Sbir(7)|6=o = sinh y/^jL/ ^7 which is 
the result for a wire with Dirichlet boundary at its ends. The reservoir acting as a phase breaker, 
the phase sensitivity is lost since the reservoir is now located in the ring. 

Neumann boundary at the end of the arm.— It is also interesting to study the spectral 
determinant for Neumann boundary at the end of the arm, Ai = 0, which describes a reflection 
of the diffusive particle (isolated network). In this case [7] 

5 N cu(7) = 2 cosh ^jb [cosh V7^Sr cu " cos 6] (16) 
where the effective perimeter now reads 

cosh ^/t"L^ u = cosh y/^yL + - sinh ^p^L tanh ^b . (17) 

In this case the limit b — > leads to the spectral determinant of an isolated ring £Neu(7)|&=o = 
2(cosh ypyL — cos 9). 

We now analyze the different behaviours in a time representation to understand their physical 
significance in the diffusion problem. To analyze the diffusion at time t, we have to consider 
7 = 1/t. In the following we keep using the notation 7 = 1/Lz, for convenience. In addition 
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to the Thouless time tl = I? needed to explore the ring, there is a second characteristic time 
Tfe = b 2 which is the typical time required to explore the arm. 

• Short times t <C tl, t& (i.e. -C L, b).— In this case the arm is explored over a distance 
smaller than the perimeter and the presence of the arm has a small influence. The precise 
nature of the boundary condition is not felt when diffusing inside the loop, since the trajectories 
encircling the loop do not have enough time to reach the end of the arm (t <C T&). This is 
reflected by the fact that L^ r = L^ u in this limit. From the expressions (|15I17[) . we obtain 
the effective parameter L c g ~ L + L !fi ln(3/2). We recover the behaviour of the isolated ring 
n t ~ ^exp-^. 

• Intermediate times tx t t\, (i.e. L <C <C b).— The particle can turn diffusively 
many times around the loop but cannot explore the arm up to its end. This is a limit of an 
infinitely long arm. Like in the previous case, the boundary condition plays no role. From 
(|15I17|) . the effective length is L g ~ \fL^ l l 4 = ^JLL^. By using (|TU|) we obtain : 

dt W n (t) e-* ~ e- n ^^ 1/4 , (18) 
which leads to the harmonic content (jBJ) 2 ■ The inverse Laplace transform gives 

»-.M-»M^x(<-=£) (w) 

where x(£) = | / °° due - "*"^" cos (-j=£u + f ) = | Re(e i f / °° due"^")) with p(u) = -u 4 + 

■u£e _1 4. #(t) is the Heaviside function. The function at the origin is x(0) = ^^/^ while it 
presents an exponential tail : 

x(0 - 4= 777^173 e " 3 (5/4)4/3 fOT f » 1 • ( 2 °) 

V6vr (£/4) i/M 
As a result, the tail of the distribution behaves like 

1 n 4 / 3 

w »® ~ T^iT^ exp - C 7W ' (21) 

where c is a coefficient. We will see in section [4.41 that for a fixed initial condition W n (x, x;t) 
presents the same exponential behaviour with a different prefactor. This expression shows that 
the winding number scales like nt oc 

• Large times t^,tt, <C t (i.e. L,b <C Lin).— In this regime the diffusive trajectories can 
explore the arm until its end and the precise nature of boundary condition matters. 

For Neumann boundary condition, the expansion of ((T7j) gives L^ n ~ \/L(L + b). It is 
interesting to point that since the arm is explored until its end, this result can be simply obtained 
by replacing L v by L+b in the result L c r ~ ^JLL^ obtained for L <C <C b. Since this effective 
perimeter does not depend on 7, it describes normal diffusion around the ring. However the 
diffusion constant is reduced due to the time spent in the arm nt — y/2t^' 2 /yfl? + bL. 

For Dirichlet boundary, from (|15jl. the effective length reads ~ L^sjLjb. This reflects a 

behaviour f °° dt e~~ /t W n (t) oc e~ n V L / b independent on 7 which originates from the absorption 
at vertex 1. The winding number reaches a finite value at infinite time : nt — yJlbjL. 

We summarize these results in tabled 

2 Note that only the exponential behaviour and the exponent of 7 in the prefactor coincide in © and 1181 
since this latter has been obtained by a uniform integration of the cooperon over the network. 
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4.2 Relation with the anomalous diffusion along the skeleton of a comb 

We discuss specifically the regime tl -C t <C t\, of an infinitely long arm. From (|19j) we have 
seen that the winding number scales with the time like 

n t oc t 1/4 . (22) 

The prefactor can be obtained from the method explained in the appendix[Bl We recall that the 
winding in a ring without arm (normal diffusion) is nt — ^/2t 1 / 2 / L. In the presence of the arm, 
the Brownian trajectories can explore the arm over long distance when turning around the ring. 
For a time scale t the effective perimeter of winding trajectories is L e g (7 = 1/i) ~ y/Lt 1 ^. From 
this simple heuristic argument we recover the subdiffusive behaviour nt ~ t 1 ^ 2 /L c g ~ i 1//4 /\/L. 
We stress that "sub diffusion" refers to the time dependence of the winding number around the 
flux, and not to the motion inside the wires which obeys normal diffusion. 

This problem is similar to the known problem of the diffusion along the skeleton of a comb : 
the diffusion in the ring is the periodisation of the diffusion along the skeleton of the comb (see 
figure [U). This problem has been studied by a different method in refs. [S] and [§] (note that 
the reference |§] corrected a wrong assumption about the nature of the distribution made in 
[H]), however the power law in front of the exponential was not given. The problem of diffusion 
along the skeleton of a comb belongs to a broader class of problems : the diffusion of a particle 
along a line with an arbitrary distribution of the waiting time r spent on each site. It was 
shown that if the distribution of the waiting time presents an algebraic tail P\{t) oc t -1- ^ with 
< fi < 1, the diffusion is subdiffusive with nt ~ W 2 [EJ- I n the case of the diffusion in the 
comb, the distribution of the trapping time by the arm is given by the first return probability 
for the one-dimensional diffusion : P\{t) oc t -3 / 2 and we recover nt ~ i 1 / 4 . 



4.3 The case of N a arms 

By using the mapping between the diffusion inside the ring and in the comb we immediately get 
the result for N a arms : the trajectory encounters N a arms for one turn, which corresponds to 
N a steps of length L/N a in the comb. Then the result for N a arms is obtained by a substitution 

n -> nN a and L -» L/N a which leads to J °° dt W n (t) ~ e -nv^L7 1/4 . For N a = 2 

the exponential dependence agrees with (JHJ). 

• N a arms attached regularly.- The result given by the previous simple argument can also be 
obtained from a calculation of <S(7) for several arms. If _/V a arms of length b are attached regularly 
around the ring we obtain 



ra=l 



cosh(-^— ) - cos(^^ ) + - smh(— ) coth ^76 



(23) 



for Dirichlet boundary conditions at the end of the arms. For iV a = 1 we recover ()14|) . (The 
spectral determinant for Neumann boundary conditions is given in chapter 5 of JJ)- F° r short 
times, t <C tl^t^, we get L c g ~ L + N a L^ ln(3/2). The term N a L v ln(3/2) has the same origin as 
the term (2/3) 2n in Q : it is related to the probability to cross the A^ a vertices of coordinence 3 
for a trajectory of winding n = 1 7 . For intermediate times tl *C t <C t^, the expansion of the 
spectral determinant shows that L e g ~ yjNaLL^, in agreement with the above discussion. For 
large times T£,T{, <C t, the effective length reads L e fj ~ L v ^/ N a L/b that reflects the absorption 
of the particle by the reservoirs after a finite time. For N a = 2, that is for the ring of figure El 
we give the expression of the effective perimeter : 

cosh ^/7-L e ft = cosh y/jL + sinh y/^L coth yp^b + - sinh 2 (^/7L/2) coth 2 yf^b . (24) 



8 



• N a arms attached at the same point. - Note that the study of the effect of several arms is easier 
if we consider the case of N a arms attached at the same point. It immediatly follows that 

/ s i n h pyb \ Na / N \ 

5(7) = 2 ^ — J j^cosh - cos 9 + -y sinh coth ^ybj (25) 

for Dirichlet conditions, which leads to cosh^/7~L c fj = cosh^yL + ^ sinh ^/7-L coth^/76. In 
the short time limit, t <C tl,T{,, we obtain : L e g ~ L + L ¥ ,ln((A r a + 2)/2), whose second term 
originates from the fact that each winding requires to cross one vertex of coordinence N a + 2 
[7j. For large time t 3> T£, the winding properties are exactly similar to the case of regularly 
attached arms. 

4.4 Fixing the starting point 

In this paragraph we study the distribution of windings with a fixed starting point xq. We choose 
the origin of the Brownian trajectory to be at the vertex for simplicity (see figureEH)- Following 
H'd.2\ we consider the new spectral determinant S^ x °\j) for generalized boundary condition with 
parameter Ao at the vertex 0. It is easy to obtain : 

S^\,) = S(,) + X Smh ^ Lsmh ^ b . (26) 

7 

By using we immediately get the cooperon at the vertex 

1 sinlu/7L , 

p c (0,0) = (27) 

2^/7 cosh y/^yLcs — cos 6 



f 

Jo 



from which 

' dt W n (0, 0; t) e-T* = J- Sl f f T L e-Vf^ . (28) 
In the limit L <C <C 6, it reads 

dtW n (0,0;i)e^ t ~^ I7i e-^ 1/4 . (29) 
The inverse Laplace transform gives 

W n (0, 0; t) = OWjfli *(t= ( 3 °) 

where ip(£) = ^ Re(e _1 4 duu 2 e~^ u ^). The function <p(u) is defined above. At the origin we 
have "0(0) = ^^/^ while i[>(Q presents an exponential tail : 

V>(£) ^ -4= (£/4) 1/3 e" 3(5/4)4/3 for £ > 1 • (31) 
V6vr 

The tail of W n (0, 0; i) then reads 

^1/3 j^4/3 

W n (0,0;t) oc^exp-c-^. (32) 

The exponential behaviour is the same in Q2UJI and (|31[) while the exponent of the prefactor 
changes because we have fixed the initial condition instead of averaging over it. 
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5 Ring connected to a network 



The above examples show that the harmonics of the AAS oscillations of a mesoscopic ring can 
change drastically due to the excursion of Brownian trajectories in the arms connected to the 
ring. An interesting question is how the harmonics are modified if the ring is connected to a 
network more complex than a wire ? An important point is to know whether the Brownian 
motion inside the network is recurrent or not (a Brownian motion is said to be recurrent if 
it comes back to its starting point with probability one after an infinite time). The fact that 
the one-dimensional random walk is recurrent is crucial to lead to the behaviour (|6I18I29|) . as 
explained in section FOl In dimension d > 2 the Brownian motion is known to be not recurrent 
whereas the case of d = 2 is marginal. Therefore we expect the dimension 2 to play a special role. 
This is one of the reasons why we consider below the case of a ring connected to a square two- 
dimensional network, as shown on the figure The formalism is introduced for an arbitrary 
network. 




Figure 4: A diffusive ring attached to (a) : an infinite square lattice, (b) : an infinite Bethe 
lattice of coordinence z = 3. In both cases the bonds have equal lengths a. 

First we consider the network without the loop. It is characterized by a matrix whose 
determinant gives the spectral determinant 5 ne t(7). 

If we now attach a ring at the vertex of the network, the new network is characterized by 
a matrix M = M nc t + SAi with 

6Map = A(9,Xo)5 M 6po (33) 

where (appendix C of [7j) 

. , . s ~- cosh ,/7-L — cos 9 , . 

AC,A.)-A. + Vr £^_ L . (34) 

Since 5M has only one nonzero element, (&M)oo> we have 

detM = [1 + A(0, A ) (M~^)oo\ det7W net . (35) 
This shows that the spectral determinant of the network with the loop is now 

5 (Ao) (7) = smh ^ L [i + A (9, Ac) (ACJoo] S n M . (36) 
V7 

Using (fTTj) we have P c (0, 0) = ^ln[l + A(9,X ) (-A4~ e t)oo]| Ao=0 that leads to the same result 
Q27I28JI as when the ring is connected to an arm. It is interesting to note that the matrix element 
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is related to the Green's function of the diffusion operator in the network, i.e. the cooperon at 
the vertex in the absence of the ring : 



P c net (0,0) = <0| -|0> = (.M-i)oo • (37) 

7 — A 

Therefore we can write the cooperon at the vertex in terms of the cooperon of the network in 
the absence of the ring : 

1 sinh yHL , s 

p c (0, 0) = — — . (38) 

2 \H cosh ypyL + \ J^pnctJ 0fi) ~ cos 9 

Now the effective perimeter is related to the matrix Ai nc t characterizing the network in the 
absence of the loop : 

1 sinh A /7~L 

cosh /rLeff = cosh a/7-L + - , (39) 

2 V / 7(" M net)00 

which is a generalization of the results (|15I17|) obtained for the ring attached to an arm. We can 
easily check that if the network is simply a wire with Dirichlet boundary at one end, we have 

(A4net)oo = l/(v / 7 cotn v / 7 6 )' leadin g to <H3>- In the llmit \fl L < 1 («- e - L < L <^)> e Q- © 
leads to 

^"\GE = \CSy- (40) 

The case of regular networks.— It is interesting to consider the specific case of a regu- 
lar network. In this case we can write : -M ne t = sin ^^ a A'" where the matrix N is N a p = 
6 a p z cosh y/^a — a a p. The matrix a a p is the connectivity matrix (see appendix EJ. We have in- 
troduced the coordinence z of the network. In this case -M~ e * is related to the (discrete) Green's 
function of the connectivity matrix G(a,@) = (N^ 1 )^, defined by ^^(Eda^ — a afl )G([i, (3) = 
6 a /3, for an "energy" E = z cosh y/^a. It follows that the effective length is related to the discrete 
Green's function of the network at the point where the ring is attached 

1 sinh^/rTL 

cosh ^L cS = cosh JyL + - sinhv ^ aG(Q)Q) • (41) 
In the limit y/^jL < 1, we have 



5.1 Ring connected to a square network 

We now apply these results to the case of a square lattice. The problem left is to estimate the 
discrete Green's function of the square network (figure 0^). We have 

i r +w d 2 k i i / i \ 

' 2 J_ n (27r) 2 2 cosh y/ja — cos k x — cos k y 27rcosh^/7a Ycoshy^yay 

where K(x) is the complete elliptic integral of first kind. We have a new characteristic time 
r a = a 2 needed to explore the bond. We now consider two limits of interest: 

• Intermediate time r^, < f < t„ (i.e. L -C <C a).- In this limit the Brownian 
trajectories starting from the ring can only explore a small portion of the 4 arms to which it is 
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connected. Eq. gives L c g ~ \/4L7 x / 4 . This result corresponds precisely to the case of a 
ring connected to N a = 4 long arms and we have / °° dt W„(0, 0; t) e"^ oc exp — n\/4L7 1 ' 4 - 

• Long time Tjj • t a t [i.e. L,a<^. L^).- In this case the Brownian trajectories encircling 
the flux n times can explore the square network over distances much larger than a. Eq. (|43j) 
leads to 

^^^J—^^r. (44) 



4n(4/^7a) ' 



The probability then reads 



[ it W n{0 , 0; «) * J£ .,,(4/^) exp -nj — ^ . (45) 

The Brownian trajectories encircling the ring can leak over long distances in the square network. 
The effective diffusion in the ring is even slowed down compare to the case of a ring connected 
to arms. The number of windings behaves like : 



n t-\J^Z lnt ( 46 ) 

(the way to obtain the precise coefficient is explained in the appendix El). An argument similar 
to the one of section II~2*1 can be developed to obtain by different means. Since the diffusion 
is recurrent in the square network, this latter acts as a trap in which the diffusive particle stays 
during a time r and eventually come back inside the ring. The distribution of the trapping time 
is given by the first return probability in a square lattice, that is known to behave at large times 
like P x (t) oc l/(rln 2 r) [Hj 3 . It can be shown from this distribution that the winding number 
scales like n t oc v / fn7 [13] . 

The result ([45 1) can be interpreted as the amplitude of the n-th harmonic of the AAS oscil- 
lations of the conductivity : 



Ao- (n) oc exp -n, / 1 . . for L, a < L v . (47) 

y ahi(4L v /a) v 

5.2 Ring connected to networks of higher dimensions 

To emphasize the role of the dimension of the network, let us consider a d-dimensional hyper- 
cubic network. The Green's function reads G(0, 0) = i j* BZ [d cosh yj^a — ^f =1 cos fcj] -1 , 
where the integration is performed over the Brillouin zone. This can be conveniently rewritten 
(5(0,0) = \ J °° dy [e~ ycosh ^ a Io(y)] d where Io(y) is the modified Bessel function, (i) In dimen- 
sion d = 1 and d = 2, the integral is dominated by the neighbourhood of k ~ (or the domain 
of large y) in the limit y^a <C 1. We can write G(0,0) ~ Jbz ffbj? W+^dh or — 
\ J°° dy j2~yi/2 e~ yia2d / 2 . In dimension d = I the integral diverges as G(0, 0) ~ l/a^/j, which 
reflects the recurrence of the Id Brownian motion and leads to y/jL e g ~ 7 

The di . 

mension d = 2 is the marginal case : the integral diverges logarithmically G(0,0) ~ ln(l/a^/7) 
which still indicates a recurrent Brownian motion and gives (|44[). (ii) In dimension d > 2 the 
integral reaches a finite value G(0, 0) = | J °° dy [e _y /o(y)] d = fid 111 the limit 7 — > 0. Note that 
Pd — 1 / (2d) for large dimensions, that coincides with the result given below for a Bethe lattice 

3 In this work the authors considered the survival probability Wabs(r, t) for a particle diffusing from r in the 
presence of an absorbing site at 0. The probability for the particle to reach the origin for the first time is 
proportional to —dtW a bs(r,t). 
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of coordinence z = 2d. The Brownian motion is not recurrent. Therefore ^/7~L e ff is independent 
on 7, which indicates that the winding number n± reaches a finite value in the infinite time limit. 

Bethe lattice.— Let us consider the case of the Bethe lattice of coordinence z (figure Efr), that 
models a network with an infinite effective dimension. The Green's function at coinciding point 

is [m 

G(0,0) = ifb(!) with F (x) = 2(z-1) = =■ (48) 

E \EJ z - 2 + ^Jz 2 - A(z - l)x 2 

We obtain the following behaviour for t S> tl, T a 



o 



(1/ IT",, (0. 0: / ) e-~ ft oc exp -n\j ^ _ ^ ^ . (49) 



This result is similar to the one obtained for a ring connected to an arm with absorption at the 
end (Dirichlet condition). In this latter case the absorption at the end of the arm is responsible 
from the fact that the particle leaves the ring after a finite time. In the Bethe lattice, since 
the diffusion is not recurrent, the diffusive particle injected in the ring eventually get lost in the 
infinite lattice. At large times the number of windings in the ring reaches the finite limit : 



2(z-l) a . 

n ^\jW^2)L- (50) 

The higher the coordinence, the smaller the time spent by the particle in the ring is. The shorter 
a, the faster the particle feels the structure of the Bethe lattice and get lost in it. 

All results are summarized in tabled 



6 Conclusion 

We have studied the winding properties inside a ring connected to a network. Our analysis was 
based on the introduction of the effective perimeter L e s (7) that probes the winding at time scale 
t ~ I/7. We have obtained this effective perimeter as a function of the matrix «Mnet describing 
the network that is connected to the ring. The analysis of the effective perimeter immediatly 
gives the nature of the winding (normal or anomalous) since the winding number scales with 
time as rit ~ yt/L e g(l/t). 

Our study of winding properties was motivated by the physics of quantum transport through 
a ring. We have emphasized the importance of taking properly into account the external wires 
connecting the ring when studying quantum transport. Experimentally, the measurement of ra- 
tio of harmonics provides a direct way to extract the phase coherence length and its temperature 
dependence (this has been used very recently on large square networks ^H])- In particular assum- 
ing an exponential behaviour exp —nL/L^,, given by ©, or the nonexponential exp —ny/2L/L ip , 
given by l(fi)l. leads to very different temperature dependences of L„(T). The figure El summarizes 
our main result : we have plotted the effective length as a function of L„. The linear behaviour 
at small L v , corresponds to © and the square root behaviour at large L v , to (jjjj). Since the 
crossover between the two regimes occurs on the range L ~ L v , it is useful to know the nature 
of this crossover. We note that the effective perimeter for the ring with N a arms can be well 
approximated by L e g ~ \f L 2 + N a L v L which interpolates between L and N a L v L (see inset 
of figure 03). Therefore the harmonics decay approximatively as exp — n^J (L/L^) 2 + NaL/L^. 
It is clear from this approximation that the crossover, that occurs for L„ ~ L/N a , can be more 
easily reached for a large number of arms. In a metallic ring, the phase coherence can reach 
several microns. Therefore, the regime L v > L/N a seems to be attainable experimentally. 
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Network 


Regime 


^cfr(7) 


winding raj? 


no 


t <C Tjj 
T L < t 


L 
L 


normal 2 exp — 4r 
normal 2t/L 2 


N a arms 

Neumann 
Dincnlet 


t < t l , n 

tl, r b < i 


L 

^/NaLL^ 

y/L(L + iV a &) 

r / AT T / U 

L^^NaL/b 


normal 2 exp — hr 
anomalous 2^/iri/{N a L) 
normal 2t/(L 2 + N a bL) 
limited 2b/{i\ a L) 


square 


r L < t < r a 

TL, T a <^t 


T 1 2tvL 


anomalous \Art/(2L) 
anomalous m (i) 


d-hypercubic 
a > I 


t l < t < r a 
TL, Ta <C t 


^2dLL v 
T 1 1 L 


anomalous \/7ri / (<££) 
limited zprfj 


Bethe 


r L < t < r a 
TL, Ta < * 


yj zLL^ 

T Jz(z-2)L 
z-1 a 


anomalous 2y/irt/(zL) 
limited ? 

2(2—2) 1/ 



Table 1: We summarize the results for the winding around the ring of perimeter L connected 
to a network. We recall that 7 = l/L^ allows to probe the various time regimes since it is the 
parameter conjugate to time : 7 ~ 1/t. There are in general three characteristic times : the time 
tl = L 2 to diffuse around the ring, the time T a = a 2 to explore one bond of the network and the 
time Tfc = b 2 to reach the boundary of the network of linear size b. Except for the case of the 
wire, we have considered infinite networks with b = 00. 

3 
2 

1 





12 3 

W L 

Figure 5: The effective perimeter L e g for the ring of figure^ with N a = 2 arms, as a function 
of the phase coherence length. L e g is given by eq. \2J$ . We have taken the limit b 3> L^. The 
dashed line is the expansion for large phase coherence length L c g = y / 2L ip L(l + + •••). 
The thin line is the linear expansion L c r = L + ln(9/4) + • • • at small L v . In the inset, we 
compare the exact result with the approximation \JL 2 + 2L ip L (dotted line). 

The phase coherence length has been introduced as an effective parameter 7 = 1 jl?^ in the 
cooperon P c (x, x') = (x\ | x' ) in order to describe phase breaking mechanisms. However, the 
effect of electron-electron interactions, that dominates at low temperature, is not well described 
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by such an effective parameter |16| . It has been shown recently that the behaviour of the AAS 
harmonics in a ring behaves in this case as exp — nT 1//2 L 3 / 2 in the limit <C L , where T is 
the temperature ^7j. However this behaviour has not been observed in the recent experiments 
|15j on large square networks, in which the condition <C L is not well fulfilled. A mechanism 
similar to the one discussed in our paper could explain the discrepancy between the theory and 
the experiment : for L v > L the excursion of the Brownian curves of finite winding in the 
surrounding part of the network increases the effective perimeter of these paths and modifies 
the dependence of the harmonics in L/L^. 



Acknowledgements 

C. T. acknowledges stimulating discussions with Alain Comtet and Jean Desbois. We are grateful 
to Satya Majumdar for his enlighting remarks and having pointed to our attention Ref. |12j . 



A Spectral determinant 

A particularly convenient way to describe the spectrum {E n } of the Laplace operator A on a 
network of one-dimensional wires is to consider the spectral determinant £(7) = det(7 — A) = 
Ylni'y+En) where 7 is the spectral parameter. If the loops of the network are pierced by magnetic 
fluxes the derivative is replaced by the covariant derivative : A — > D 2 with = -j=- — iA(x). 

Let us consider a network of B wires connected at V vertices. These latters are labelled with 
greek indices a, /?,... First we introduce the connectivity matrix a Q/ g characterizing the topology 
of the network : a a p = 1 is a are (3 connected by a bond (we denote (a/3) the bond and l a p 
its length). a a p = otherwise. If we define a scalar function ip(x) on the network, boundary 
conditions at the vertices must be specified. At the vertex a we choose (i) continuity of the 
function, that is all the components ip a /3(x) of the function along the wires af3 issuing from a 
tend to the same limit as the coordinate reaches the vertex, (ii) a a /3~D x ijj a p(a) = X a ^(a) 
where the connectivity matrix in the sum constrains it to run over the neighbouring vertices of 
a. Therefore the sum runs over all wires issuing from the vertex. The real parameter X a allows 
to describe several boundary condition: A Q = 00 forces the function to vanish, ip{a) = 0, and 
corresponds to Dirichlet boundary condition that describes the connection to a reservoir that 
absorbs particles. \ a = corresponds to Neumann condition, which ensures conservation of the 
probability current and describes an internal vertex. 

It was shown in J] that the spectral determinant is : 

S(7 ) = n **A<(7) (51) 

(a/3) ^ 7 

where the product runs over all bonds of the network. M is a V x V matrix defined as : 

K + > ] a ai i coth( A /7^ At ) - a a p v — - , (52) 
~ / smh(^/77 a/3 ) 

where the connectivity matrix constrains the sum to run over all vertices /i connected to a. 
We have included the magnetic fluxes {9 a p is the flux along the wire). For a more detailed 
presentation, see jj]. 



B Time dependence of the winding number 

We explain how to compute efficiently the time dependence of the winding number around a 
ring connected to a network. For this purpose we compute (n 2 )t, where the average is taken for 
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all close trajectories starting from a specified point and coming back to it after a time t 



2 _ E n w 2 W n (a,a;t) 
( >' " £ n ^(*,*;t) (53) 

((n)t = follows from the symmetry P c \g = P c \_g). 

To go further we consider the case where the initial condition x is the vertex where the 
ring is attached to the network. Then the cooperon has the structure (|27|). Let us define the 
quantity 

1 d 



From © it follows that f2 m (7) = J °° dt e~ 7 * ^ n n m W n (0, 0; t) which is related to the m-th 
moment of the winding number. We have 

2 _ r- 1 [n 2 ( 7 )] (55) 

{n }t ~ £-i [o (7)] (55) 

where £ _1 [•■•] designates the inverse Laplace transform. From (|27|) we can obtain the two 
general expressions : 

^o(7 = ~ri^~f \2 ( 56 

4 V7 ( smh y^) 2 

rw ^ 1 sinh^L 

2 (7 = j- . 4 ■ (57) 

8^7 ( sinh 



Example 1 : the isolated ring. - The effective perimeter is equal to the perimeter in this case : 
L cS = L. We have O ( 7 ) = ^ coth(^7L/2) and f2 2 ( 7 ) = ^ coth(yfyL/2)/ smh 2 (^fL/2). 
We consider two regimes : 

• Short times t <C tl '■ We have ^0(7) — jyf and ^2(7) — ^ e ~^ L which gives after Laplace 

transform {n 2 )t — 2exp— This result was obvious from the expression (|13[). 

• Long times t tl : We obtain in this case ^0(7) — ^£ and ^2(7) — which leads to 
(n 2 ) t ~ 2t/L 2 . This is the normal diffusion. 

Example 2 : the ring connected to a square network. - We consider the long time limit t ^> Ti,T a 
to demonstrate (|46[). The effective perimeter is given by (|44|). We find ^0(7) — 4if m (16/a 2 7) 
and ^2(7) — g^|; ln 2 (16/a 2 7). Using the fact that in the limit 7 — > we have £ _1 [In I/7] ~ 1/t 
and £ _1 [in 2 I/7] ~ 21n(i)/i, we eventually find (n 2 ) t ~ ^ Int. 
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